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1. INTRODUCTION 


The OGO-VI satellite, which was launched on June 5, 1969 carried 
a complement of twenty-six experiments. One of those instruments, the 
F-02 package, was a cylindrical Langmuir probe experiment whose primary 
purpose was to measure ionospheric electron temperatures and densities. 
This report briefly describes the F-02 experiment itself, outlines the 
computer programs developed at The University of Michigan to analyze 
the raw data, and gives a summary of the scientific information obtained, 
so far, with the aid of this experiment. 
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2. EXPERIMENT DESCRIPTION 


Cylindrical Langmuir probes have been used successfully for nearly 
a decade to measure electron temperatures and densities in the ionosphere 
(e.g. Brace, et al. , 1963; Taylor, et al . , 1963; Nagy, et al. , 1963; 
Spencer, et al. , 1965; Brace, e_t al . , 1971). The theory of cylindrical 
Langmuir probes was first published nealy half a century ago by Mott— 
Smith and Langmuir (1926). More recently, this theory was extended to 
moving probes by Kanal (1964) . The rocket- and satellite-borne instru- 
mentation associated with these probes has been discussed extensively 
in the open literature (see Brace, et al., 1971). 

The two cylindrical probes used on OGO-VI were mounted on the OPEP 
at right angles to each other as indicated in Figure 1. Each probe as- 
sembly consisted of a collector 9 in. long and 0.022 in. in diameter, 
and a concentric guard 9 in. long and 0.065 in. in diameter. The guard 
and the collector were made of stainless steel, and teflon was used as 
an insulator between them. The potentials of both the guard and the 
collector were driven together, but only the current to the collector 
was monitored. The purpose of the guard was twofold: 1) to keep the 

collector physically away from the OPEP surfaces, and 2) to help main- 
tain a uniform field around the collector. A large portion of the 
OPEP surface was covered with a fine mesh gold-plated grid in order 
to insure a uniform sheath about the OPEP and to increase the "current 
dump" area of the spacecraft. 

The electronics package was housed in the main body of the space- 
craft. It consisted of a two-range logarithmic current detector and 
ramp voltage generator system (System I), a four-range linear current 
detector and ramp voltage generator system (System II), and a common 
control system shown in Figure 2. Tables I and II give many of the 
details of the characteristics of Systems I and II respectively, while 
Table III gives the command status indications. 

The current detector in System I is a two-range differential one 
which has the following transfer functions: 

%ut - 2 lo %0 (3160 X p + » 
e out ' 2 lo8 10 < 31 ' 6 + » 
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where e^^ and I are in volts and microamperes, respectively. The de- 
tector ranges alternate with each ramp voltage sweep. The ramp voltage 
is a linear 6 V sweep with a period of 9.2 sec for normal operation and 
2.3 sec in the fast format mode of operation. 

A spacecraft will acquire some variable negative potential with re- 
spect to the plasma. The probes have to be swept both positive and neg- 
ative with respect to the plasma and therefore it is desirable to be able 
to adjust the starting voltage of the ramp voltage sweep. In System I 
a special circuit which adjusts the ramp voltage bias is used to achieve 
this. At the beginning of each sweep the bias is automatically adjusted 
until the detector output is zero. The ramp voltage then sweeps from 
-2.5 V to +3.5 V with respect to this zero current (floating) potential. 
Thus the magnitude of the sweep is always 6 V but the starting point is 
allowed to vary from -4.5 V to +7.5 V, depending on the spacecraft po- 
tential. 

The output of System I consists of three 144 msec (36 msec in fast 
format) bits followed by a voltage which is proportional to the logarithm 
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of the collector current. The three 144 msec bits carry the following in- 
formation: 1) ramp voltage retrace, 2) range, probe, and calibration in- 

formation, and 3) the ramp start voltage. The system was normally oper- 
ated with a 9.2 sec ramp voltage period with measurements alternating 
sequentially between probes 1 and 2. The system could be switched to 
fast format (2.3 sec period), and/or to measurements with one of the 
probes only. 

An in-flight calibration took place approximately every ten minutes. 
During the calibration sequence a known resistor was attached to the in- 
put of the current detector to simulate full scale detector current. 

During this calibration period a complete ramp voltage sweep was also 
placed on the data output line to allow the ramp voltage to be monitored 
directly. 

System II is quite similar to System I in operation but much simpler. 
System II utilizes a four-range differential current detector and produces 
as its output a voltage proportional to the collected probe current. The 
appropriate transfer functions are the following: 
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where e QUt and Ip are in volts and microamperes, respectively. The de- 
tector changes ranges after each ramp voltage period which is 9.2 sec 
under normal operation and 2.3 sec during fast format operation. 

This system does not have an automatically self-adjusting ramp 
starting voltage like that of System I; in this system the ramp voltage 
sweeps linearly from -3.5 V to 2.5 V with respect to a bias voltage 
level (one of the six bias levels, namely 0 V, 2 V, 4 V, 6 V, 8 V, 10 V, was 
selected by impulse command during the operation of System II). The for- 
mat of the data output from tills system is identical to that of System I: 
three information bits followed by probe current for each sweep. The 
bits carry the same information as indicated for System I. The calibration 
period for this system also occurs approximately every 10 minutes, and con- 
sists of zero current checks followed by detector current outputs due to 
appropriate resistors being placed across the input terminals in order to 
simulate full scale currents. During the final portion of the calibration 
sequence the ramp voltage is monitored directly. 

Data from the sub com provided all the needed information for the 
impulse and power command status of both experiments. The details of 
the subcom monitor are shown in Table III. 
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Figure 1. OGO-F-02 probe mounting configuration. 
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TABLE I 


SYSTEM I CHARACTERISTICS 


Sensitivities (Full Scale) Range 1 0.1 x 10 6 A 


Range 2 10 x 10 6 A 


Power 28 Vdc 2.24 W 
(21.5 to 50 Vdc) 

AV Normal Format -2.5 V to +3.5 V + Bias* 

— a PP 

= *682 V/sec 
dt 


Fast Format 

= 2.728 V/sec 
at 


-2.5 V to +3.5 V + Bias* 


Sweep Characteristics 



Detector Characteristics 


Range 1 


0.1 yA full scale 

e = 2 log ' {3160 V/yA|l (yA)| + l} 
out 10 p 


*Bias is V for zero net probe current ±100 mV; it is searched at each 
app 

retrace time. 
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TABLE I (Continued) 


Range 2 

10 yA full scale 

e = 2 log.. 31.6 V/yA| I (yA) | + 1} 
out 10 p 

Output Characteristics 

Data Output 

z = 10 3 0 . 
o 

Data Limited at 4.95 V 


Format (Time scale exaggerated) 


Bit 


9.216/2.304 sec 


Bit 1 SYNC BIT, ALWAYS 5.1 V 
Bit 2 RANGE, PROBE, CALIBRATE BIT 


VOLTAGE • 


FUNCTION 



0.5 

MEASURE 

PROBE 2 

RANGE 

2 

1.0 

II 

" 1 

II 

2 

1.5 

II 

" 2 

II 

1 

2.0 

II 

" 1 

It 

1 

2.5 

CALIBRATE 

" 2 

II 

2 

3.0 

11 

" 1 

II 

2 

3.5 

II 

" 2 

II 

1 

4.0 

II 

" 1 

II 

1 


5.0V -f 


Bit 1 


Bit 2 




144/36 msec 


0.0V 


-t= 
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TABLE I (Concluded) 


Bit 3 


SWEEP VOLTAGE START BIT 
FLIGHT UNIT: 

SWEEP VOLTAGE START = 2.666 (Bit 3) 


- 4.909 



TABLE II 


SYSTEM II CHARACTERISTICS 


Sensitivities (Full Scale) Range 

1 

0.01 

yA 

Range 

2 

0.07 

yA 

Range 

3 

0.5 

yA 

Range 

4 

3.5 

yA 


Power 28 Vdc at 1.88 W 
(21.5 to 50 Vdc) 


AV -3.5 V to +2.5 V + Bias* 

— ■ *EB- 

dV 

-j- = .682 V/sec 
dt 

*Bias is 0V, 2V, 4V, 6V, or 10V depending on command status 


Sweep Characteristics 



Bit Period 


:i< 



(432 msec) 
Active Period 


(9.216 


.432 sec) 


Detector Characteristics 


Range 1 


0.01 ]iA full scale 

e - 1.54 V + 350 V/yA(I (yA)) 
out p 
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TABLE II (Continued) 


yA full scale 

= 1.54 V + 50.0 V/yA(I (yA)) 

P 

yA full scale 

= 1.54 V + 7.01 V/yA(I (yA)) 

p 


yA full scale 

= 1.54 V + 1.00 V/yA(I o (yA)) 

Output Characterisitcs 

Data Output 

z = 10 3 n 
o 

Data Limited at 4.95 V Negative Clamp at -0.6 V 

Format 



Range 2 

0.07 

e , 
out 

Range 3 

0.5 

e , 
out 

Range 4 

3.5 

e 

out 


0 
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TABLE II (Concluded) 


Bit 1 SYNC BIT ALWAYS 5.1 V 
Bit 2 RANGE, PROBE, CALIBRATE BIT 


VOLTAGE 

0.500 

0.750 

1.000 

1.250 

1.500 

1.750 

2.000 

2.250 

2.500 

2.750 

3.000 

3.250 

3.500 

3.750 

4.000 

4.250 


PANGE 4 
" 3 

" 2 
" 1 
" 4 

" 3 

" 2 
" 1 
" 4 

" 3 

" 2 
" 1 
" A 

" 3 

" 2 
n x 


PROBE 2 



MEASURE 


CALIBRATE 


COMMON MODE CK 


M 


Bit 3 SWEEP VOLTAGE START BIT 
FLIGHT UNIT: 

SWEEP VOLTAGE START = 2.174 (Bit 3) - 3.620 
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TABLE III 


SYSTEMS I AND II 

COMMAND STATUS 

INDICATIONS 




SUBCOMMUTATOR 

OUTPUT LEVELS 





VOLTAGE 


FUNCTION 


CMD 

STATUS 

0.500 

SYSTEM I 

RESET 


1 

1 1 

0.750 

1 ! 

FAST MODE 

PROBE, 2 ONLY 

1 

1 0 

1.000 

fl 

it ii 

ii ii 

1 

o“r 

1.250 

It 

it ti 

ALTERNATE 

1 

0 0 

1.500 

tv 

RESET 


0 

i i 

1.750 

It 

NORMAL MODE 

PROBE 2 ONLY 

0 

1 0 

2.000 

T! 

ii ti 

ii j it 

0 

0-1- 

2.250 

II 

ii n 

ALTERNATE 

0 

0 0 

2.500 

SYSTEM II 

RESET 


1 

1 1 

2.750 

It 

AV BIAS=10V 


1 

1 0 

3.000 

II 

" " 8V 


1 

0 1 

3.250 

fl 

" " 6V 


1 

0 0 

3.500 

II 

RESET 


0 

1 1 

3.750 

II 

AV BIAS= 4V 


0 

1 0 

4.000 

II 

ii ii 2V 


0 

0 1 

4.250 

M 

" " ov 


0 

0 0 
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3. DATA ANALYSIS 


This section outlines the computer program which was developed 
and used to take the raw data from the decom and orbit tapes provided 
by the project and 1) reduce the ionospheric electron temperature and 
density as well as the spacecraft potential results from the F-02 ex- 
periment, and 2) to display these results along with other supporting 
information (e.g., angle of attack, latitude, longitude, etc.)- Before 
discussing the program and the numerical techniques, it is appropriate 
to outline briefly those aspects of the theory of moving cylindrical . 
Langmuir probes which were used by the data analysis programs to ob- 
tain atmospheric electron temperatures and densities. 

3.1. Probe Characteristics 

3.1.1. Infinite Probes 


The classical theory describing the current collection (volt-ampere) 
characteristics for stationary probes of simple geometry is well-known 
(Mott-Smith and Langmuir, 1926). When an infinitely long cylindrical 
probe whose radius is less than the characteristic Debye length is im- 
mersed in a plasma, the current collected is — t— _ 


I 


A Ne 
c 



exp (— V) 


v < v 
P 




/V + e V erfc(/v)| - A N e_/-~ . ~i_(l-+-V) 

/it 


-v— v- 


( 2 ) 


where V = I e (v-v ) /kT I 
P 

I = current due to the species under consideration 
e = electron charge 
v = applied potential 


14 



= space potential (plasma zero reference) 

N = density 

A = collector area 
c 

k = Boltzmann's constant 
T = temperature. 

This theory has been extended (Kanal, 1964) to moving probes in 
Maxwellian plasmas: 


I = A Ne 
c 



exp [- (V+k 2 ) ] 


CO ^ 

7 ■ (k//v) I ( 2 k/v) 

0 (n!) 2 2k n 


v <v 

P 


(3) 


where k = XsinQ 

X = u//2kT/m 
u = spacecraft velocity 
0 = angle of attack 

I = modified Bessel function of first kind 
n 

n 

(k "nT ~ F (n + 3/2,V)J n (2k/v) v >v p (4) 


I = 



V-k* 


l 

0 


where J is the n tk Bessel function and T (n + 3/2, V) is the incomplete 
n 

gamma function. 

The 0G0 spacecraft orbital velocity (approximately 8. km/sec) and a 
mean thermal speed for electrons of approximately 300 km/sec result in a 
speed ratio of about 0.025. For these values of X Equations 3 and 4 differ 

from Equations 1 and 2 by negligibly small amounts ( 0.1%). Thus the probes— -- 

may be considered stationary with respect to electrons. 

However, the heavier ions are moving considerably more slowly such 
that the speed ratio is about 4, so that the random motions of these par- 
ticles are less important than the drift component. Indeed, if one ignores 
the thermal motion, assuming a cold ion model (the limiting model as X -*■ °°) , 
the following equation describes the collected current (v > v ) 
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(5) 


I 


A’Ne 

c 




e(v-v ) 

El 

u 


where A' = cross-sectional area of the probe 
c 

U = %mu 2 . 


When Equation 5 is compared to Equation 4 (X = 4), once again, the. 
differences are negligibly small (<0.1%), so that the thermal motions 
of the ions may be ignored. 

Thus, a complete model describing current collected by a long thin 
probe moving through a plasma from both ions and electrons is as follows: 


I = 


I = 


= A Ne 
c 

/kT 
1 e 

e(v-v p ) 1 

1 / 1<Ti 

V 2 V exp 

kT j 

! \/ 2m . it 

v X 

= A Ne 1 
C 1 

1 /kT 

1 / e 

e (v-v ) \ 
P 1 

/ kT. 
_ / i 

[\/ 2m e tr exp 

kT 1 

\j 2nu 7T 


e(v -v) 

A'Ne \/l 4 — v < v ; X > 4 

c V U p — 


1 + 


'±sl"r\ 

kT I ! 


v < v ; X = 0- 
P 


/kT ? / e (v-v ) 

I = A Ne »/-=— •— \ /l + — + ion 
c V 2mir yjf V kT 


V > V 


where the ion component when v > v p is ignored, because the v = v p ion 
current is a factor of 50 smaller than the electron component and is 
decreasing as v increases. 


3.1.2. Finite Probes 
3. 1.2.1. Retarding Region 

The F-02 experiment had two probes mounted on the orbital plane 
experiment package so that one probe was always at 90 0, to the velocity 
vector and one at 0° to the velocity vector. The probe mounted at 90° 
to the velocity vector was found to conform well to the above-outlined 
theory when v < v . However, because of a well-known end effect, the 


(6a) 
~(6b ) - 
(7) 
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ion current collected by the other probe, aligned with the velocity 
vector, did not behave at all in accordance with Equation 4. The 
essence of the explanation for this phenomenon concerns the influence 
of the spherical sheath surrounding the tip of the finite (~9 in.) 
probe, through which an ion must pass before being collected. 

Although this end effect has been treated in the literature and is 
understood qualitatively, there is no useful quantitative description of 
this component of the collected current which can be utilized in a data 
reduction scheme. Therefore, the procedure used for F-02 data from this 
probe was to assume a linear model which adequately describes the current: • 


I = a, + a 0 v + A Ne 
12 c \ 


kT 


2m it 
e 


exp 


f e(v-v ) 
l + kT 


( 8 ) 


Examples of the use of this model and model 6a are presented in Figure 3. 
3. 1.2. 2. Accelerating Region 


In the accelerating region for electrons (v > v , so that electrons 
are attracted to the probe), neither the 0° probe or the 90° probe conformed 
to the theory. The proposed explanation for this effect once again concerns 
the "finiteness" of the probe. 

Unpublished work by M. Kanal has suggested that a semi-infinite 


cylindrical probe, i.e., 
described by 

a probe with 

one and only 

one end. 

would be 

I = A Ne 
c 

/ kT / v 

e^erfc (/v) ' 

i+v 

(9) 

v 2m7T [w * 

2 

J 

+ 2 J 


The first term in Equation^ is one-half the current which would be 
collected by an infinite cylinder; the second term is one-half the cur- 
rent which would be collected by a sphere of equivalent surface area. 
Figure 4 presents the theoretical normalized currents from a sphere, 
an infinite cylinder. Equation 9, and the result of a "best fit" power 

Q l 

law model (1 + V) , where a = 0.85. It is clear that the semi- inf inite 
cylinder model describes a probe characteristic considerably better than 
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the model for an infinite cylinder. 

The results of computer calculations of accelerated current col- 
lected by a finite cylinder due to streams of monoenergetic charged 
particles indicated that, when the potential distribution about the 
probe was described by a simple power law, the current collected was 
"reasonably linear" except near space potential. Inasmuch as a Max- 
wellian distribution may be viewed as a superposition of a number of 
monoenergetic beams, it is not unreasonable to assume that the current 
collected due to Maxwellian particles passing through the spherical 
sheath about the tip of the probe might also be reasonably linear. 
Because of this, and for lack of a better alternative, the following 
model, with modified linear term, was used in the accelerating region 
for the finite probes aboard OGO-VI. 


An example of the correspondence between the data and this model is 
shown in Figure 5. Although the derivation of this model is empirical 
and somewhat arbitrary, it is far superior to the infinite cylinder, at 
least for the F-02 experiment. 

3.2. General Program Flow 

Figure 6 is a block diagram describing the program flow. The F-02 
data reduction program consists of a main routine (MR) and several sub- 
routines, MR has two major functions: 1) it acts as an interface be- 
tween the person requesting data reduction and the machine, and 2) it 
directs the program flow to the currect subroutines, MR first checks 
to see of the requested outputs are consistent with each other. If 
not, it makes an appropriate change to whatever data is to be output. 

The routine then loads the subroutines needed to perform the requested 
data reduction. All output tapes are then put in their appropriate 
starting positions, and the data request to be reduced is read from 
cards. 



/v e erfc(/V) 
tt 2 


1+pV 

2 


( 10 ) 
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MR is now ready to begin its second function. Subroutine LABEL is 
called to read the decom tape label of the file to be processed. The 
label is decoded and printed and all subroutine setup entries are called. 

If part of a file is to be reduced, rather than the whole file, its limits 
are read and the tape is moved to the starting time. Subroutine SWEEP is 
called to gather the next voltage sweep and resulting data curve, and the 
analog bits are decoded. Subroutine CALBRT is called to determine if the 
curve is part of a calibration sequence. If so, new calibration constants 
are computed and EXPPRT is called to print the ones that changed. If the 
curve is not a calibration curve, subroutine CONVRT is called to convert 

the input data to volts and amperes utilizing the standard relationship 

Y = mx + b, and this subroutine in turn calls the data reduction subroutine 
BILL which reduces the raw probe data to temperature. When subroutine 
BILL returns to CONVRT, CONVRT prints the reduced data along with selected 
orbital parameters (Figure 7) . 

When all of these subroutines return to MR the subroutine INPRT is 
called to print any experiment parameters which changed for this par- 
ticular curve, and also to print appropriate comment about its disposition 
(Figure 8). At this point the program returns to get the next data sweep 
and this continues until the requested stop time is passed or an end-of- 
file condition on the input tape is detected, error conditions notwith- 
standing. At this point the program will reduce another file on the 
current tape, or on another tape, or stop, depending on the availability 
of another data request and its contents. 

In addition to the tasks of producing numbers, the data analysis sub- 
routines also provide for the optional display of data curves and resulting 
fits in four different forms. The experimenter has his choice of viewing 
his data and results as 

1) a listing of all the raw data and computed results for each curve, 

2) a printer plot of the raw data and the "residuals" from a fit, 

3) a Calcomp plot of the raw data and the residuals, 

4) a storage tube CRT display which can be filmed, 

5) any combination of the above. 
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Figure 3. Graphical program output, exhibiting retarding 
region data and model equations, for both probes. 
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NORMALIZED CURRENT 



Figure 4. Theoretical normalized currents vs. normalized 
potential for a sphere, an infinite cylinder, Equation 9, 
and power law model (1 + V) a , a = 0.85. 
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COLLECTED CURRENT (jjA) 


2.0 


1.5 



3.5 4.0 4.5 5.0 5.5 


APPLIED VOLTAGE 

Figure 5. Comparison of Equation 10 with data. 
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Figure 7. Sample of reduced data output. 
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Figure 8. Spacecraft parameters. 



3.3. Data Fitting Technique 


Once raw data has been gathered from decom tapes and converted to 
voltages and current, subroutines BILL arid either TMPRTR, TMPBI, or DENSTY 
are called to obtain values for the atmospheric parameters, density, tem- 
perature, and space potential. The technique applied can be called "Least 
Squares Fitting," "Maximum Likelihood Estimation," or "Nonlinear Regression." 
In any case the data collected during a voltage sweep are viewed as a random 
sample from a probability density whose mean is a function of applied vol- 
tage which we refer to as a "model". Thus: 


Y 'ri N(f(a,v i ),a 2 ) 


or 

Y = f(a,v ± ) + e ± 'ri N (O,0 2 ) 

where Y . is a random variable representing current collected at voltage v^ 
v^ is the voltage applied to the probe 

e. is a random variable representing the "error" or noise associated 
1 with the observation of current Y 

a 2 is the variance of the noise distribution 

^N(u,a 2 ) means normally distributed with mean u and variance 0 2 

a is a vector of parameters to the model 

f(a,v.) is the modeling . function, or a form which the data would follow 
1 as a function of voltage in the absence of noise, according to 
the best available theory (Equations 6a, 8, and 10). 

Esti ma tes of the parameters to the model are computed by assuming 
that the proper values are those which would make the observed "random 
sample" the most likely one. In the event that the "errors" are normally 
distributed this condition is equivalent to a least squared error crite- 
rion, which is useful even when errors are not normally distributed. Thus 
statistics for the atmospheric parameters are calculated by finding the 
vector a such that 
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min 


0 = 1 ( Y , - f (a,v_.)) 2 = 
i=l 


or 


„ n 3f(a,v ) 

(Y.-fCa,v.))— ^-0 
3 1=1 3 


j = 1, • • -,m 


where m is the number of unknown parameters to the model and n is the 
number of observations of current and voltage in the random sample.. 

Because of the nature of the models for current collected by the 
F-02 probe, the above system of equations is intrinsically nonlinear. 

Hence an iterative numerical technique must be used. The one chosen here 
is a modified Newton's method for systems of nonlinear equations which is 
implemented by subroutine STNLSE. It is equivalent to the "Linearization 
of the nonlinear model on an iterative basis. Although more specialized 
techniques have been applied to probe data (Theis, 1968), 'this approach 
is used because it is more general (three different models are used), it 
converges slightly faster, is easily modified, and does -not entail a sig- 
nificant amount of extra computation. 

In order to prevent an inadequate model in the accelerating region 
from biasing the results obtained from a good model in the retarding 
region, the two regions were reduced separately. The algorithm used to 
locate these two regions takes advantage -of the fact that the difference 
between space potential (v ) and the floating potential (or zero current 
potential, which is the most salient feature of a data curve) when ex- 
pressed in units of normalized potential (i.e., multiplied by e/kT^) is 
theoretically a weak function of temperature with a value of approximately 
2. Thus the floating potential is determined and the voltage equivalent 
of 1.2 units of normalized potential at a temperature which is an average 
of the three previous temperatures is added. Then all data from the 
beginning of the voltage sweep to this voltage are considered to be the 
retarding region: 

{ (Y ,v )|v < v. < v-, + 1.2 kT/e} = retarding region 

i l 1 s — x — f 1 
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where I = instantaneous current at voltage v^ 

= first voltage in sweep 
v^ = instantaneous voltage at index i 
v c = voltage such that Y (v c ) = 0 

T = average of three previous temperatures 
lc = Boltzmann's constant 
e = charge on an electron 

Similarly, the accelerated region is defined as: 

3 ?kT 

{(Y.,v.)| v,, H < v. <v } = retarding region 

i i 1 fl e i e a a 

where v^ = last voltage in sweep. 


3.4. Error Estimates ” - - 

Since the least squares estimates used to determine values for the 
atmospheric parameters in the above models are themselves random variables, 
they have complicated probability distributions whose variances can be ... - 
estimated. These variances are quite useful as error estimates because they 
can be used to establish confidence intervals for temperature estimates at 
any desired level of risk. However, since the models are nonlinear, the 
level of risk can only be approximated and is accurate only to the extent 
that the model can be considered linear within the confidence interval. 
Nevertheless, approximate error estimates are almost as useful as accurate 
ones, and are certainly more useful than none at all. Therefore, subroutine 
STDERR was used on all temperature fits to calculate the variance of the 
model parameter from which temperature is computed. From this variance, 
the width of an approximate 80% confidence interval in units of degrees 
Kelvin was calculated. This number is returned along with computed tem- 
peratures to the main program and provides a convenient indicator of data 
quality. 
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4. RESULTS 


A preliminary version of the data analysis program was operational 
on the Goddard Space Flight Center computer at the time of launch, and 
therefore some near-real-time data analysis was carried out during the 
first ten days of the satellite's lifetime. 

The first presentation of data from the F-02 experiment was given 
informally at the December 1969 AGU Meeting. A special session of the 
1970 National AGU Meeting held in Washington D.C. was devoted to pre- 
sentations of early results from the OGO-VI satellite, and a paper on 
F-02 results from June 1969 was presented at that session (Samir, Nagy, 
and Brace, 1970). At the April 1971 National AGU Meeting another 
special session was organized for the presentation of OGO-VI data, this 
time for data related to the geomagnetic storm of November 1969. Re- 
sults from the F-02 Langmuir probe experiment were presented in two of the 
papers (Reber, et al., 1971; Hanson, et al, 1971). 

O 

The winter predawn enhancement of the 6300 A oxygen red line 
emission and the electron temperature has received a great deal of 
attention in the last few years. The F-02 experiment along with the F-03 
ion trap experiment have provided excellent information on the electron 
temperature increase and on the integrated intensity of conjugate photo- 
electrons. These results have been presented at the April 1971 National 
AGU Meeting by Sanatani, et al. (1971) . 

The most recent theory of red arc formation is that given by Cornwall, 
Coroniti, and Thorne (1970), who suggest that the ring current protons 
generate ion cyclotron waves in the plasmapause region, and that a fraction 
of this wave energy is transferred to plasmaspheric electrons via Landau 
damping. This energy, via thermal conduction, is transported down to the 
ionosphere causing the red line emission. The 0G0-VI satellite carried 
a complement of experiments capable of measuring atmospheric parameters 
that are relevant and necessary for a complete understanding of the pro- 
cesses involved in the formation of a red arc. Because of a number of 
factors (e.g. spacecraft potential problems, infrequent red arc occurrences, 
etc.), only one red arc has been studied in detail, the results of which 
have been presented (Nagy, et_ al. , 1972a), and published (Nagy, et al. , 1972b) 
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The OGO-VI experimental results were in excellent agreement wi th all as- 
pects of the theory with the possible exception of the proposed wave- 
particle interaction mechanism. 

Discrepancies between electron temperatures measured by Langmuir 
probes and radar backscatter techniques have been noted in recent years 
(e.g. Carlson and Sayers, 1970). A concentrated effort was made to 
study this discrepancy with the aid of the F-02 probe experiment. All 
the available comparison points between electron temperatures measured 
by the F-02 probe and the various backscatter stations, as well as the 
ion temperature comparisons between the F-03 ion trap and backscatter 
data, have been collected, The results of these comparisons have been 
presented (Nagy, et al. , 1971), and published (McClure, et' al., 1972). 

The average discrepancy between the probe and radar-deduced electron 
temperatures was found to be about 12%, or less than most of the previous 
comparisons. There were insufficient statistics available to subdivide 
the data according to important theoretical parameters (e.g. probe diam- 
eter — Debye length ratio, temperature, etc.) and therefore no clear 
conclusions regarding the reason for the discrepancy were possible. 

Electron and ion temperature results from the F-02 and F-03 ex- 
periments, respectively, have shown very significant temperature minima 
in the nighttime equatorial region at altitudes above about 500 km. A 
great deal of effort was made to ascertain that these unexpected tem- 
perature minima were real and not a result of instrumental effects. 

When the integrity of the data was established it became clear that 
these low electron and ion temperatures are indications of the hitherto 
unexpected importance of adiabatic cooling of the plasma in the equa- 
torial region. The results of these observations and the related theory 
are in the process of being published (Hanson, Nagy, and Moffett, 1972), 
and further work on this subject is in progress. Simultaneous electron 
and ion temperature observations by F-02 and F-03 have indicated also 
that ion temperatures may exceed electron temperatures at night (Hanson, 
Nagy, and Moffett, 1972), due to the higher electron thermal conductivity 
as pointed out earlier by Nagy, Bauer, and Fontheim (1968). 



5. FUTURE ACTIVITIES 


In order to "minimize cost and maximize science" the data re- 
duction has proceeded on a very selective basis. Data from some 
selected orbits are being reduced at the moment to study the equa- 
torial adiabatic cooling phenomenon further. It is also anticipated 
that studies of the angle of attack behavior of the probes will be 
carried out in the near future. The Langmuir probe experiments on 
the ESR0-1A satellite have indicated that the electron temperature 
is not isotropic in the nighttime midlatitude ionosphere at higher 
altitudes. An attempt will be made to see whether this phenomenon 
can be studied with data from the .F-02 experiment. 
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